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ABSTRACT 

Among the newly discovered features of multiple stellar populations in Globular Clus- 
ters, the cluster NGC 1851 harbours a double subgiant branch, that can be explained 
in terms of two stellar generations, only slightly differing in age, the younger one hav- 
ing an increased total C+N+O abundance. Thanks to this difference in the chemistry, 
O ■ a fit can be made to the subgiant branches, roughly consistent with the C+N+O 

S—{ ' abundance variations already discovered two decades ago, and confirmed by recent 

| spectroscopic data. We compute theoretical isochrones for the main sequence turnoff, 

by adopting four chemical mixtures for the opacities and nuclear reaction rates. The 
standard mixture has Z=10~ 3 and [a/Fe]=0.4, the others have C+N+O respectively 
equal to 2, 3 and 5 times the standard mixture, according to the element abundance 

■ distribution described in the text. We compare tracks and isochrones, and show how 
' the results depend on the total CNO abundance. We notice that different initial CNO 
, abundances between two clusters, otherwise similar in metallicity and age, may lead 

t ' to differences in the turnoff morphology that can be easily attributed to an age dif- 

ference. We simulate the main sequence and subgiant branch data for NGC 1851 and 
show that an increase of C+N+O by a factor ~3 best reproduces the shift between the 
' subgiant branches. According to spectroscopic data by Yong et al., the C+N+O abun- 

. dance in this cluster appears correlated with the abundance of s-process elements, Na 

and Al, and this makes massive AGBs the best progenitors of the C+N+O enriched 
population. We compare the main sequence width in the color 1x1^336 w-nrFsiw with 
models, and find that the maximum helium abundance compatible with the data is 

■ Y~0.29. We consider the result in the framework of the formation of the second stellar 

generation in globular clusters, for the bulk of which we estimate a helium abundance 
of Y ^ 0.26. The precise value depends on which are the AGB masses from which the 
C+N+O enriched matter originates, and on the amount of dilution with the pristine 
gas. 

Key words: globular clusters: general; globular clusters:individual: NGC 1851; 
stars : abundances 



1 INTRODUCTION 

The observations of Globular Cluster (GC) stars still need 
to be interpreted in a fully consistent frame. Neverthe- 
less, a general consensus is emerging on the fact that most 
GCs can not be considered any longer simple stellar pop- 
ulations, and that "self-enrichment" is a common feature 
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among them. The first suspicions originated from the well 
known "chemical anomalies", already noted in the seven- 
ties (such as the Na-0 and Mg-Al anticorrelations) . Re- 
cently observed to b e present at the turn o ff (TO) and among 
the s ubgiants (e.g. iGratton et~ai1 1200 ll ; iBrilev et ail |2002| . 
120041 ) . they must be attributed to some process of self- 
enrichment occurring at the first stages of the cluster life. 
A decisive feature indicating the presence of more than 
one population in GCs has been the observation of mul- 
tiple, well separated seque nces. Multiple main sequences 
have been found in u Cen l|Bedin et alj |2004| ; iPiotto et alj 
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Figure 1. Variation with temperature of the opacity of the mix- 
tures used in the present investigation. At the top of the figures, 
the value of R=density/temperature 3 and of the hydrogen mass 
fraction X are labelled. Upper panel: opacity values for the stan- 
dard (solid track), a— enhanced case, and for the mixtures with 
CNO abundances increased by a factor 3 (dotted) and 5 (dashed), 
corresponding to a total metallicity Z=0. 00235 and Z=0. 00350. 
For the two latter cases we plot the values for three total metal- 
licities Z=l, 2, and 4xl0 — 3 (from left to right). Lower panel: 
opacities in the low-temperature regime, for the mixtures with 
CNO enhancements of 2 (solid), 3 (dotted), and 5 (dashed). 



20051) and NGC 2808 (jD'Antona et alJ 120051 : IPiotto et all 
20071 ). Besides, multip le subgiant branch es have been ob- 



served in NGC 1851 (jMilone et ail 120081) an d other clus- 
ters, among which NGC 6388 (|Piottol 120091 ) . The former 
phenomenon can only be interprete d in terms of popula- 
tions with d ifferent helium content l)D'Antona et al.ll2002l ; 
lNorrisll2004l ). T he latter may be discussed in terms of a dif- 
ference in age l|Milone et al.l l2008) . or in total CNO abun- 



dance (jCassisi et al.ll200sl ). and is the subject of the present 
work. In both cases, we are left with the difficulty of finding 
a consistent origin for these populations. We can speculate 
that there was a first epoch of star formation, that gave ori- 
gin to the "normal" (first generation, hereinafter FG) stars, 
with CNO and other abundances similar to Population II 
field stars of the same metallicity. Afterwards, there must 
have been some other epoch of star formation (second gener- 
ation, hereinafter SG), including material heavily processed 
through the hot CNO cycle in the progenitors, belonging to 
the FG, but not enriched in the heavy elements expected in 
supernova ejecta. This material either comes entirely from 
the stars belonging to the first stellar generation, or it is a 
mixture of processed gas and pristine matter of the initial 
star forming cloud. 

The SG in mo st clusters is a high fraction of the to- 
tal nu mber of stars (|Carretta et alJlboOfit iD'Antona fc Carol 
2008). In order to have enough CNO processed material 
available, it is necessary that "self-enrichment" is a result 
of pollution from an initial stellar population much larger 
than the stellar content of today's GCs. This is generally 
attributed either to the dynamic loss of a great fraction 
of the clusters' FG s tars in the early evolutionary phases 
jP'Ercole et al.ll2008l ), or to the formation of the GC within 
a much larger stellar environment, such as a dwarf galaxy, 
where the polluting matter is supplied by the surrounding 
stars (|Bekki fc Norrisll2006l ). 

One of the important constraints for the progenitors 
was usually considered to be that their matter must have 
been processed through the hot CNO cycle, and not, or 
only marginally, through the helium burning phases, since 
the sum of CNO elements is the same within observa- 
tional errors in the "normal" and in the anomalous stars 
l|Cohen fc Melende dl2005l ; llvans et alJll999h . If the polluting 
matter is identified with the envelopes of asy mptotic giant 
branch (AGB) stars (| Ventura et alj|200ll . 12003 ). these AGBs 
must be very massive, so that their evolution is only scarcely 
affected by the third dredge up, that acts to increase the 
surface abundance of the primary carbon formed in the he- 
lium intershell during the thermal pu lses, and then partiall y 
mixed into the external envelope (e.g. llben fc Renzinilll983l ). 
Otherwhise, pollution must come from a totally different 
kind of objects, such as the envelopes of fast rotating mas- 
sive s tars during the core H-burning phase l|Mevnet et al.l 
120061 ; iDecressin et al.l [20071 ) . and in this case we do not ex- 
pect any C+N+O enhancement at all, as mentioned before. 

Recently, the constancy of C+N+O has been challenged 
by the discovery that the clus ter NGC 1851 disp lays a dou- 
ble subgiant branch (SGB) dMilone et all [2008). No such 
splitting is present in the main sequence of NGC 1851, in- 
dicating a normal — or close to normal (see Sect. 4) — he- 
liu m content fo r both SGBs. This has been interpreted 
by I Milone et alJ (|2008l ) as due to a ~lGyr age difference 
between the two SGB populations. This age difference is 
much larger than ge nerally believed possible for a GC (e.g. 
iD'Ercole et alj|2008l). An alte rnative interpretation has been 
given by ICassisi et al.l (|2008l ): they were able to reproduce 
the splitting by assuming that stars in the faint SGB (fSGB) 
have a larger C+N+O abundance and similar age of the 
bright SGB (bSGB). At the same time, lYong fc Grundahl 
(2008) have shown that this cluster harbors a star-to-star 
abundance variation in the s-process elements Zr and La 
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and that the abundances of these elements were correlated 
with the abundances of Na and Al, and anticorrelated with 
O. Furthermore, within the small sample of 8 giants exam- 
ined, there was a hint that the abu ndances of the s-p rocess 
elements is bimodal. Most recently. lYong et al.l (|2009h have 
determined the C, N, and O abundances in four red gi- 
ants, selected to span the range of Na, Al and s-process 
abundances, and found that indeed the sum of C+N+O 
exhibits a range of 0.6dex, a factor 4, and that the light 
elements abundances, s-process and total C+N+O abun- 
dances are correlated. These observations support indeed the 
AGB stars as progenitors of t he stars popu l ating the fSGB 
in NGC 1851, as suggested bv lCassisi et all (120081 ). Signs of 
a milder C+N+O enrichment in other clusters have been no- 
ticed bv lCarretta et all (|2005l ) in NGC 6397 and NGC 6752, 
but NGC 18 51 is to date the most prominent example (see 
|Piottoll2009l . for further examples). As we have noticed, in 
most GCs the (possible) AGB progenitors of the SG must 
be so massive that the chemistry of the ejecta is scarcely 
affected by the third dredge up, with only a minor increase 
in the total C+N+O abundance. But the effects of the third 
dredge up increase with time: i ts efficiency increases when 
the evolving mass decreases fe.g. lVentura fc D'Antonall2008l . 
Table 1). So the observability of a C+N+O increase depends 
on how long the p hase of star formation of the SG stars lasts. 
lYong et al.l (|2009h propose that in NGC 1851 this phase lasts 
long enough that the effect of the third dredge up is evident. 

In this work, we consider the evolution of low mass stars 
evolving at GC ages, with the metallicity of NGC 1851, but 
considering th e effect of an increas ed C+N+O abundance on 
the evolution. ICassisi et al.l ((2008) limited the analysis to a 
chemistry in which the C+N+O content was approximately 
doubled with respect to the standard composition, while we 
consider C+N+O abundances increased by factors 2, 3 and 
5. Selection of the individual abundances for C, N and O is 
explained in Section 2. Inclusion of their effect is considered 
both in the low temperature and high temperature opacities 
and in the nuclear reaction network. In Section 3 we present 
the results, and discuss the time evolution of selected evo- 
lutionary tracks, the total mass lost along the RGB, and 
the relative location of the isochrones in the color magni- 
tude diagram (CMD). In Section 4 we present simulations, 
in the CMD, of the main sequence and SGB for NGC 1851, 
under several hypotheses concerning the different composi- 
tions necessary to explain the separate SGBs. In Section 5 
we discuss the results. 



2 INPUT PHYSICS AND MODEL 
COMPUTATION 

All the evolutions presented in this work have been calcu- 
lated by means of the ATON code for stellar evolution, with 
the numerical structure described in details in Ventura et al. 
(1998). Tables of the equation of state are generated in the 
(gas) pressure-temperature plane, according to the OPAL 
EOS of the Livermore group (see OPAL w ebpage, last up- 
date in February 2006, iRogers et al.|[l996l ). replaced in the 
pressure ionization regime by the EOS by Saumon, Chabrier 
& Van Horn (1995), and extended to the high-density, high- 
temperature domain according to the treatment by Stoltz- 
mann & Blocker (2000). 



2.1 Standard and non-standard opacities 

For the "standard" models we adopt the latest opac- 
ities by Ferguson et al. (2005) at temperatures lower 
than 10000 K and the OPAL opacities in the version 
documented by Iglesias & Rogers (1996). The mixture 
adopted is alpha-enhanced, with [a/Fe] = 0.4 (Grevesse 
& Sauval 1998). Electron conduction opacities were taken 
from the WEB site of Potekhin (see the web page 
http://www.ioffe.rssi.ru/astro/conduct/ dated 2006) and 
correspond to the Potekhin et alj (1 19991 ) treatment. The 
electron opacities are harmonically added to the radiative 
opacities. For this project we select three different mix- 
tures of elements, having the C, N and O abundances var- 
ied with respect to the standard mixture^. For the abun- 
dances of C, N and O we adopt the values in Table 2 in 
I Ventura fc DAntonal (|2008l ). corresponding to the yields of 
the models with 5, 4.5 and 4Mq and metallicity Z=10 -3 . 
These abundances are reported in Table 10 As we see, in 
the 5 Mq ejecta Oxygen is depleted by 0.46 dex (the original 
[O/Fe] is +0.4 in the standard mixture), Carbon is enhanced 
by 0.13 dex and Nitrogen is enhanced by 1.7 dex. The total 
C+N+O is a factor 2.1 larger than in the standard mixture. 
The 4.5 Mq total CNO content is 3.1 times the standard 
one, and the 4M© one is a factor 4.9 larger. We remarked 
that t he 4 red giants of NGC 1851 examined by lYong et alj 
(2009) show a maximum variation of ~4 in C+N+O. In- 
dividually, carbon and oxygen vary up to about a factor 3 
and nitrogen by a factor 7. With our present choices, we 
are then enhancing the nitrogen increase with respect to 
these few observational data. More observations and new 
models may be useful in future analyses. The total "metal- 
licity" Z in mass fraction for the three CNO-enhanced mix- 
tures is listed in Table 1. The radiative opacities for these 
CNO enhanced mixtures have been computed on purpose 
for this work. For temperatures above 10000 K we utilized 
the online computations from the OPAL group found at 
http://physci.llnl.gov/Research/OPAL/opal.html; for lower 
temperatures we computed opacities with the Ferguson et 
al. (2005) code for the Z mixtures listed in table 1. Figure [T] 
shows the differences among the opacities of the various mix- 
tures in the high- and low-temperature regime. In the up- 
per panel we see that the largest differences are found in the 
ioniziation region of the CNO elements, around logT ~ 6. 
The lower panel of Fig. 1 shows the low temperature opaci- 
ties for the three enhanced computations. Above log T ~3.5 
few differences are noted as the opacity is dominated by hy- 
drogen. However, at lower temperatures the opacity signifi- 
cantly deviates for each run from the standard. This devia- 

1 We do not inc l ude v ariation in the sodium abundance, as done 
bv lCassisi et al" I [120081 ) , considering that its total abundance is in 
any case very low, and its contribution both to energy generation 
and opacity is negligeable. 

2 The choice of abundances has been motivated by the naive hy- 
pothesis that th e SG may directly be formed by the ejecta of the 
FG, and that the Vent ura fc D'Antonal (120081 ) yields represent the 
SG composition. In fact, in Sect. 5, we will reinterpret the results 
in the light of a dilution model. Notice however that direct forma- 
tion from the FG ejecta does not necessarily imply an anomalous 
IMF peaked at the masses whose yields are adopted, if we allow 
for "self-enrichment" from a much wider cluster environment, as 
discussed in the Introduction. 
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Table 1. Chemistry of the models 



Name 


Z 


total CNO 


[C/Fe] 


[N/Fc] 


[O/Fe] 


Y 


M/Mq(AGB) 


age (Myr) 


CNOxl (standard) 


0.00100 


1 


0.00 


0.00 


0.40 


0.240 






CNOx2 


0.00185 


2.1 


0.13 


1.70 


-0.06 


0.324 


5.0 


103 


CNOx3 


0.00235 


3.1 


0.12 


1.89 


0.19 


0.310 


4.5 


128 


CNOx5 


0.00350 


1.9 


0.14 


2.02 


0.44 


0.281 


4.0 


166 



tion is due to more and more O available (due to the CNO 
enhancements) for the formation of molecular water at these 
temperatures. The low temperature opacities do not affect 
the structure of the models we are considering, while at the 
larger temperatures of the stellar interiors some differences 
appear. In addition, the CNO total abundance affects the 
time evolution as soon as the stars begin evolving at the 
turnoff, where the CN cycle becomes important. 

2.2 Color —T fj conversions 

Our aim is to compare t he th eoretical predictions with 
the data by iMilone et all (|2008l ) for NGC 1851, obtained 
in the HST Advanced Camera for Surveys (ACS) filters 
F814W and F606W. Therefore we convert our L,T off val- 
ues into these magnitudes by using the transformations 
for ACS bands b y Bedin et a l. (200 5), based on the solar 
scaled models by Cassisi et al.l l|2004t ). We note that the a— 



enhanced versus solar-scaled transformations, for the metal- 



licity and colors we are using, do not differ for T e 



5000K. 



On the other hand, we are also dealing with CNO en- 
hanced mixtures, and in principle the use of the same trans- 
formations could not b e adequate. However, we refer to 
iPietrinferni et all ([2009) , who widely discuss this issue, and 
conclude that, in broadband filters not bluer than the John- 
son V band, the effect of this inconsistency should be very 
small, since at least the color - T e g transformations are 
largely independent on the metals and their distribution (see 
lAlonso et al]|l996l . [l999l : ICassisi et"ai1l2004r ). This statement 
luckily applies to our case, as we are dealing with red and 
near infrared bands. Possible differences in the conversions 
due to the enhanced CNO may of course change the quanti- 
tative conclusions of our analysis, so that we raise the prob- 
lem of the computation of proper color-T e ff conversions for 
CNO enhanced mixtures. 



2.3 Synthetic CMD Simulations 

For a better comparison of the models with the data, we use 
simulations for the main sequence and subgiant branch (es). 
These are obtained by extracting random each star location 
along a specified isochrone, according to a choice of the ini- 
tial mass function (IMF). We assume a power law IMF with 
exponent -1.5 (where Salpeter's is -2.3). The exact choice of 
the IMF is inconsequential, as the comparison of the SGBs 
involves a very small mass interval (see Sect. 4). Gaussian 
errors consistent with the observational errors on colors and 
magnitudes close to the turnoff are attached to each extrac- 
tion. For the double stellar population of NGC 1851, we 
use a standard isochrone of CNOxl for 55% of the sample, 
while the rest is extracted from the isochrone having the 
same age and enhanced CNO. These percentages are based 
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Figure 2. Luminosity versus time evolution for the masses 0.8, 
0.85 and 0.9 Mq (from right to left) as a function of the CNO 
enhancement. The full (black) lines are the CNOxl evolution. The 
CNOx2, CNOx3 and CNOx5 tracks are located progressively at 
smaller luminosity. 



on the observational analysis by IMilone et afl {2009), and 
are confirmed in the analysis of the subgiant branches made 
in Sect. 4, describing the comparison between the simula- 
tions and data. 



3 MODEL RESULTS 

3.1 The isochrone location: dating GCs with the 
same metallicity and different CNO 

Figure [2] shows the luminosity versus time evolution of 
masses 0.8, 0.85 and 0.9 Mq for different CNO enhancement. 
The CNOxl standard tracks have the fastest evolution, be- 
cause the CNO enhanced tracks lie at progressively lower 
main sequence luminosity. In order to understand better the 
different roles of the opacities and nuclear reaction rates, 
we show in Figure [3] the HR diagram location of a track of 
0.85 Mq, having the standard CNOxl opacities, but CNOx3 
abundances (in particular the CNOx3 abundances are used 
in the nuclear network). We see that the main sequence evo- 
lution is quite similar to that of the standard CNOxl track, 
as burning occurs mainly through the proton-proton chain. 
The main sequence of the track is just slightly more lumi- 
nous than for the standard CNOxl track, and the evolution- 
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Figure 3. Comparison between the HR diagram of the standard 
CNOxl track of 0.85 M Q (full line) and of the CNOx3 track (dot- 
ted line), with the evolution of the same mass, obtained by assum- 
ing the opacities of the standard track, but the CNO abundances 
of the CNOx3 track (dashed line). 
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Figure 5. Isochrones in the plane Mp606W vs. m,pgo6W — 
m F8i4W- The standard Y=0.24 CNOxl t=12Gyr isochrone is 
the upper line. The lower lines are for the same age and different 
helium content. The SGB luminosity is only slightly altered for 
Y=0.32. 
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Figure 4. Isochrones in the plane Mp606W vs - the color 
^F606W — A^F814W- The top group of lines refers to an age of 
9Gyr, the lower set are for an age of 13Gyr. The isochrone lumi- 
nosity decreases as a function of the CNO enhancement. 



ary times are consequently barely shorter. The turnoff, how- 
ever, occurs at a smaller luminosity, as in that case the CNO 
cycle dominates. Therefore, the lowering of the turnoff lumi- 
nosity is mainly due to the effect of the abundances on the 
final phases main sequence burning, while the lenghtening 
of the evolutionary times is mainly due to the smaller main 
sequence luminosity of the models with enhanced CNO. 



Globally, the isochrones show a monotonic shift in lumi- 
nosity and color with increasing CNO. Figure [4] shows that 
the subgiant location for the same age is fainter for higher 
CNO, so that these models can in principle reproduce the 
splitting of the SGBs in NGC 1851, if the two sequences 
have different CNO content. Figure [5] shows the compari- 
son between the standard SGB for Y=0.24 and the SGBs of 
the models CNOx3 with different helium contents (Y=0.24, 
0.26, 0.28 and 0.32). We see that the SGB shift does not 
depend on the helium content, at least for this metallicity, 
while the main sequence location and the red giants become 
bluer when Y increases (see also Sect. 5). 

Considering Fig. [4] we observe that any age determi- 
nation for a GC can not ignore what is the total C+N+O 
abundance in its stars. While this is a well known theoretical 
aspec t of evolution (e.g. lSimoda fc Ibenll970l ; lBazzano et al.l 
1982), it is worth recalling it in this context, now that we 
have interesting evidence of its possible role. In fact, the 
interp retation of NGC 1851 data shows this (|Cassisi et al.l 
2008), and we also h ave other cases where it may be of rel- 
evance (|Piottol|2009l ). The total CNO abundance in the gi- 
ants of M4 (NGC 6121) is larger than the total CNO of 
those stars, in NGC 1851, that can be considered "nor- 
mal" for what concerns the abundances of CNO, s-proces s 
elements, Na and Al (see Figure 3 in lYong et al.l 12009) . 
Therefore, should we compare these two clusters, their rel- 
ative age determination should include the dependence on 
the CNO abundance. In this context, it is worth recalling 
the lar ge spread in N abun dance in all the clusters sur- 
veyed l|Grundahl et al.l 1 1999h . together with the dichotomy 
between CN strong and CN weak stars in so many GCs. 
The giants in NGC 6752 present a star-to- star abundance 
variation in N of 1.95dex (|Yong et al.ll2008l ) . and a similar 
spread is present also in the main sequence (jPasquini et al.l 
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Figure 6. The top panel shows the mass at the helium flash, 
Mfl as h as a function of the total stellar mass, for the three labelled 
helium contents and for the standard tracks. The dotted lines 
show the same for the CNOx3 tracks. Although for a fixed mass 
there is a dependence of Mfl as h on Y and CNO, the bottom panel 
shows that the mass lost along the RGB (<5M) practically does not 
depend on the helium content at a fixed epoch. On the contrary, 
at a given age, <5M is ~ 0.01 Mq larger for the CNO enhanced 
isochrones. We plot in the figure only the CNOxl and the CNOx3 
cases, the other two cases are practically superimposed to the 
CNOx3 lines. 

l200Sl : ICarretta et al.ll2005l ). Variations in C abundance may 
compensate, but some spread in luminosity at the tu rn-off 
of this cluster should be present (|Bazzano et alj|l982h . 

3.2 The mass lost along the red giant branch 

All models are evolved by assuming a mass loss rate follow- 
ing Reimers' (1977) prescription Mr — 4- 10 _13 r;_R-^ where 



the parameter rjR has been fixed to 0.3 for all the computa- 
tions. This allows us to find the dependence of the mass lost 
in the giant phase on the helium and CNO content of the 
tracks. Of course the results are of some significance only 
if Reimers' rate is actually a good description of the mass 
loss. In a relative sense, the results may be similar for other 
mass loss rates, if they depend explicitly only on the stellar 
parameters luminosity and gravity, and do not depend, e.g., 
on the CNO content. Some results are shown in Fig. [U In 
the top panel we see that mass at the helium flash, for a 
fixed initial mass, increases with Y and decreases with the 
CNO enhancement. We plot the CNOxl and CNOx3 case, 
but the results for the other CNO rich models are similar. 
However, what matters is the mass lost at fixed age, shown 
in the bottom panel of Fi g. El This actually do es not depend 
at all on Y (see also iD'Antona fc Caloil 120081 '). while, for a 
fixed age, all the CNO enhanced models predict ~0.01 Mq 
more mass loss, as shown for the case of the CNOx3 mod- 
els. These indications are important to predict the horizontal 
branch (HB) lo cation of sta r s with different CNO and/or he- 
lium content. Cassi si et al. I (|200Sf ) have shown that the T e ff 
location of HB models is affected by their enhanced CNO, 
due to the different relative efficiency of the shell-hydrogen- 
and core-helium- burning (se e their Figure 1, but see also 
ICastellani fc Tornambel (| 19771 ) ). 



4 COMPARISONS 

Figure [7] shows the HR diagram of NGC 1851 based on 5 
images of 350s in F606W and 5 images of 350s in F814W 
from the ACS survey by lAnderson et al.l l|200gft . Thes ex- 
posure are saturated at the basis of the RGB, so the pho- 
tometry of the RGB comes from two low exposure images 
(one 20s image in F606W and one 20s image in F814W) 
and has larger errorfl The data in Fig[7]are compared with 
the isochrones. Distance and reddening are chosen to fit 
isochrones of ages respectively 9 and 12Gyr in the two pan- 
els. We see that the giant branch location agrees better with 
the 12Gyr isochrones, but it is difficult to attribute to this 
feature a predictive value for the age, as the giant branch 
location depends on the adopted convection model. 

This comparison immediately shows that the fSGB in 
the cluster lies in between the isochrones corresponding to 
CNOx2 and CNOx3. This finding is strengthened when we 
make use of the simulations built up as described in Sect. 2.3. 
Figure [8] compares the observed HR diagram with simula- 
tions for an age of 12 Gyr. The bSGB is reproduced with a 
standard CNOxl population, while the fSGB is reproduced 
assuming a second population with CNOx2 (top), CNOx3 
(middle) and CNOx5 (bottom panel). It is evident, looking 
at the bottom panel of Fig. [8] that the CNOx5 isochrones 
are not adequate to describe the splitting of the SGBs: such 
a CNO enhancement would produce a splitting much larger 
than observed. The CNOx2 isochrones produce a very small 



3 Notice that the RGB of this cluster is broadened -or even 
split— if we observe it in t he U, B, or othe r filters in which the 
CN— bands are present (e.g. iLee et al. 2009), but the reason why 
the RGB is so broad in these observations is simply due to the 
photometric error. 
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Figure 7. Comparison of isochroncs of ages 9 Gyr (top panel) 
and 12 Gyr (bottom panel) with the HR diagram of NGC 1851. 
We apply to the data a distance modulus and dereddening in 
order to fit the "standard" top isochrone to the bSGB. The 
dashed line (black) is the CNOx3 isochrone of Y=0.26. Dis- 
tance moduli and reddenings respectively: (mp606W _ ^F606w) = 
15.72, <5color=0.05mag for the top panel (9Gyr isochrones) and 
15.44mag, 0.02mag for the bottom panel (12Gyr isochrones). The 
continuous (red) lines, starting from the top to the bottom one, 
are the isochrones for Y=0.24 and respctively CNOxl, CNOx2, 
CNOx3 and CNOx5. 



splitting, the CNOx3 isochrones are the closest to the ob- 
servations. 

For a quantitative comparison of observed and simu- 
lat ed HR diagrams we go through a procedure already used 
bv lMilone et alJ (|2009l ) . Fig.|9]illustrates this three-step pro- 
cedure for the observed ACS/WFC data. We selected by 
hand two points on the fSGB (Pij,P2j) and two points on 
the bSGB (Pi,i,,P 2 ,s) with the aim of delimiting the SGB 
region where the split is most evident. These points define 
the two lines in panel (a), and only stars contained in the 
region between these lines were used in the following anal- 
ysis. In panel (b) we have transformed the CMD linearly 
into a reference frame where: the origin corresponds to Pi,&; 



m aDS ,-M M ,= 15.44 

[! ( m eo6w- m BH») = 0oa 
t=12Gyr 
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Figure 8. Comparison of NGC 1851 data (blue dots) with simula- 
tions (green) for an age of 12 Gyr. Distance modulus and assumed 
reddening are labelled. The bSGB is simulated with a population 
having normal (CNOxl) CNO, the fSGB corresponds to a popu- 
lation with CNOx2 (top), CNOx3 (middle) and CNOx5 (bottom 
panel). The corresponding isochrones of 12 Gyr are also shown. 
The helium abundance is assumed to be Y=0.24 for both popu- 
lations. 
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Pij has unit abscissa, and both P^^b and Pa,/ have unit 
ordinate. For convenience, in the following, we indicate as 
'abscissa' and 'ordinate' the abscissa and the ordinate of this 
reference frame. The dashed green line is the fiducial of the 
bSGB. We drew it by marking several points on the bSGB, 
and interpolating a line through them by means of a spline 
fit. In panel (c) we have calculated the difference between 
the 'abscissa' of each star and the 'abscissa' of the fiducial 
line (A 'abscissa'). 

The histograms in Fig. [TO] are the normalized distribu- 
tions in A 'abscissa' for stars in four A 'ordinate' intervals. 
In the three panels of this figure we compare the distribution 
of the observed data (grey histograms) and the simulations 
(black histograms). The bSGB corresponds to the popula- 
tion with normal CNO, the fainter SGB corresponds to the 
CNO x 2 (left), CNO x 3 (middle), and CNO x 5 (right). 

If we use slightly helium enriched isochrones (up to 
Y~0.28) for the simulation of the fSGB, the results do not 
change for what concerns the SGB splitting (see Fig. 4). 
Notice however that the HB observations pose a constraint 
on the maximum helium abundance of the SG in this clus- 
ter: ISalaris et alj (|2008l ) were not able to reproduce the HB 
morphology attributing to the whole blue HB side (proba- 
bly corresponding to the same population of the fSGB) with 
models having helium content as high as Y=0.28 (and for 
larger Y the situation is worse). From the width of the MS, 
however, we can derive an independent upper limit to the 
helium abundance of the SG stars, by comparing it with 
models having different Y. 

Figure [11] shows such a comparison for eight magnitude 
bins of the MS in the color TOF336W — «iF8i4W. From the mean 
of the eight Y extimates, we derive an upper limit for the 
He spread in NGC 1851. We find AY=0.043 ± 0.003. If the 
FG has Y=0.24, the maximum possible helium abundance 
value is Y=0.290±0.003. 



5 DISCUSSION 

The computed models and the comparison with the data 
show that the CNO ehnancement required to explain the 
double SGB in NGC 1851 is of about a facto r three. This 
confir ms and puts on a more quantitative basis ICassisi et all 
(2008) results. (In the following discussion, anyway, keep 
in mind that the quantitative conclusions of our analy- 
sis are based on the use of the same color-T e ff relations 
both for CNO normal and CNO enhanced mixtures). So, 
in this cluster, we need CNO enhancement to fit the HR 
diagr am features. If we add t he sp ectroscopic informa tion 
from lYong fc Grundahll (|2008h and lYong et al l (|2009h , on 
the presence of s-process enhancement and of a global CNO 
spread among the cluster giants, we may conclude that mas- 
sive AGB polluters are the only reasonable candidates to 
form the SG stars in the fSGB in NGC 1851. 

It is possible to better quantify this statement by look- 
ing at Table 1. We choose our compositions, for the ele- 
mental abundances and opacity comp utation, from the pure 
eject a of our AGB models of Z=10" 3 (| Ventura fc D'Antonal 
2008). In the table we see that a factor three in CNO is 
reached in the matter expelled by the 4.5 Mq stars. If, how- 
ever, the matter forming the SG stars comes from pure ejecta 
only, these stars should have a helium content of Y~0.31. 



Indeed such a large helium content is to be excluded, based 
1) on the main sequence, as it does not show the width 
required by such a large Y spread (see, Figure [S] and the 
discussion in Sect. 4); and 2) on the HB morphology, both 
because the HB lacks very hot stars, and because the blue 
side of the RR Lyr reg ion is not particularly overluminous 
(e.g. ISalaris et al]|2008h . 

A possible solution is that the AGB ejecta have been 
diluted with pristine matter while forming the SG stars. 
If, e.g., the matter comes preferentially from the 4Mq, the 
ejecta would have 5 times the standard CNO, as listed in 
Table 1. Dilution with 50% of pristine matter gives the re- 
quired CNO enhancement by a factor 3. In this case, the 
starting helium abundance in the ejecta was Y=0.28. Dilut- 
ing it with 50% of matter at Y=0.24, the helium content of 
the SG stars have Y=0.26, a value not in contradiction with 
what required by the HB morphology. 

It remains to be understood why the SG took such a 
long time to be formed (according to Table 1, 166Myr is 
the age at which the 4Mq has evolved and ejected its CNO 
enriched matter into the cluster), while in most other clus- 
ters it seems to be formed much earlier, within at most 100 
Myr. In fact, Table 1 shows that the masses contributing to 
the SG must be larger than ~5Mq, if the total CNO con- 
tent has to be kept with i n a factor two of the FG value 
l|Cohen fc Melendej 120051 ; llvans et all Il999h . In addition, 
the very high h elium subpop ulations found in a few very 
massive clusters (|Piottoll20Qgf ) require that, in these clusters, 
part of the SG has been formed by the ejecta of the most 
mass ive super-AGB stars dPumo et al . 2008; D'Ercole et al.l 
l200St) . evolving at ages - 30 - 40Myr. 

The necessity of a C+N+O enhancement in the SG of 
NG C 1851 is further est ablished with our analysis, following 
the lCassisi et all l|2008h discu ssion, and has at l east an ini- 
tial observational basis in the lYong et a l. (2009) data. The 
above detailed formation scenario for this cluster, however, 
relies on the yields from the models bv lVentura fc D'Antonal 
(2008) that have been employed in this discussion. Different 
models, e.g. the set bv lKarakas fc Lattanziol (|2007t ) . provide 
stronger effect of the third dredge up at larger initial masses, 
mainly because these models adopt of a less efficient convec- 
tion d escripti on (see the discussion in Ventura fc D'Antonal 
1200581 ). The iKarakas fc Lattanziol (|2007t) yields do not 
favour the AGB interpretation of observational data for 
most clusters, namely the quasi-consta ncy of C+N+O an d 
the shape of the anticorrelation O-Na (Fe nner et alj|2003 ). 
A possible escape is that the whole range of polluter masses 
for the SG is shifted to the super-AGB range, but, to our 
knowledge, no such mo dels are available by now. H owever, 
yields close to those bv lKarakas fc Lattanziol (|2007f ) would 
change our conclusions on the timescale of formation of an 
SG, just because the stronger effect of the third dredge up 
— limited, anyway, to a factor ~3 increase in C+N+O — 
would be obtained at larger evolving masses (~ 6Mq) and 
much shorter evolutionary times. 



6 CONCLUSIONS 

We computed the evolution of standard a-enhanced low 
mass stellar models of metallicity Z=10 -3 , and of models 
with the same [Fe/H], in which the total CNO abundance is 



NGC1851 ACS/WFC data 
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Figure 9. This figure illustrates the procedure for the comparison of the observed data and the simulation. Panel (a) is a zoom of the 
observed ACS/WFC HR diagram. The two lines delimit the portion of the HR diagram where the split is most evident. Only stars 
from this region are used to measure the population ratio. In Panel (b) we have linearly transformed the reference frame of Panel (a). 
The green dashed line is the fiducial of the region around the bSGB. In Panel (c) we plotted stars between the two lines but after the 
subtraction of the 'abscissa'. 
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OBSERVED DATA SIMULATION 

Figure 10. The histograms show the A'abscissa' distribution for stars in four A'ordinate' bins. We compare the distribution of the 
observed data (grey histograms) and the simulations (black histograms). 
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Figure 11. Panel a shows the mp336W vs. rrtp336w — "1F814W CMD already presented in Milone et al. (2008) with the fiducial line 
overimposed. Panel b shows the same CMD, after subtracting from each star the color of the fiducial line appropriate for its mp336\v 
magnitude. The histograms of the color distribution in eight magnitude bins are plotted in panel c. In panel d we show the helium content 
Y as a function of the color distribution in each bin. The three black triangles indicate the MS spread predicted by teoretical models 
for a SG with CNOx3 and Y=0.26, 0.28, and 0.32, and are used to derive the straight line by means of a least squares fit. The open 
triangles corresponds to the case where the SG has CNOxl and Y=0.24. Red circles indicate the Y value corresponding to the observed 
MS color spread. 



larger by factors 2, 3 and 5, according to the elemental dis- 
tribution expected from the intermediate mass AGB ejecta 
that can be progenitors of SG stars in Globular Clusters. Our 
aim is to derive quantitative information about the C+N+O 
difference required to explain the splitting of the subgiant 
branch in the GC NGC 1851. Comparison of the models with 
the data quantifies the necessity of C+N+O enhancement 
in the fSGB. The amount of total CNO required is a factor 
of about three larger than the total CNO adopted for the 
bSGB. We incidentally warn that a difference in total CNO 
between clusters having the same metallicity may simulate 
a non negligible age difference. We try to interpret the data 
on the basis of the massive AGB scenario for the formation 
of the fSGB, showing the presence of an SG. The ejecta of 
the AGBs that formed the SG must have a total CNO abun- 
dance increased by a factor three at least. Though, if these 
pure ejecta have formed the SG, its helium abundance would 



be as high as Y=0.31, incompatible with both the MS color 
width and the morphology of the HB. We propose that the 
CNO enrichment in the ejecta was a factor at least 5 larger 
than the initial CNO, and that the ejecta have been diluted 
with pristine gas by 50%. In this case, the abundance of 
helium in the SG comes down to a much more reasonable 
value of Y=0.26 or less, but we must explain why the SG 
formation was delayed to a total age of more than 150Myr in 
this cluster and not in others. Finally, we remember that the 
quantitative conclusions of our analysis are based on the use 
of the same color-T e ff relations both for CNO normal and 
CNO enhanced mixtures. Alt hough this seems reason able 
for the red bands we are using l|Pietrinferni et aLll2009l '). the 
possibility that appropriate color-T e ff relations, when avail- 
able, may modify the CNO enahncement required to fit the 
two SGBs in this cluster must be kept in mind. 
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